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Abstract. Eigenfunctions of the Askey- Wilson second order q-difference operator for 
< q < 1 and \q\ = 1 are constructed as formal matrix coefficients of the principal series 
representation of the quantized universal enveloping algebra U q (sl(2, C)). The eigenfunc- 
tions are in integral form and may be viewed as analogues of Euler's integral representation 
for Gauss' hypergeometric series. We show that for < q < 1 the resulting eigenfunction 
can be rewritten as a very-well-poised 8</>7-series, and reduces for special parameter values 
to a natural elliptic analogue of the cosine kernel. 
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1. Introduction 

The aim of this paper is to simplify the quantum group construction of explicit eigen- 
functions of the second order Askey- Wilson g-difference operator, and to extend the results 
to the interesting and less well studied \q\ = 1 case. 

The approach is based on the known fact from [TI], and [0] that the second order 
Askey- Wilson difference operator arises as radial part of the quantum Casimir element 
of U q (sl(2,C)) when the radial part is computed with respect to Koornwinder's twisted 
primitive elements. Using this result, we construct nonpolynomial eigenfunctions of the 
Askey- Wilson second order difference operator as matrix coefficients of the principal series 
representation of U q (sl(2, C)). The two cases < q < 1 and \q\ = 1 will be treated seper- 
ately. The theory for < q < 1 is related to the noncompact quantum group U q (su(l, 1)), 
while for \q\ = 1 it is related to the noncompact quantum group U q (sl(2,M.)). 

This approach was considered in |8] for < q < 1 using an explicit realization of the 
principal series representation on Z 2 (Z). The resulting eigenfunction then appears as a 
non-symmetric Poisson type kernel involving nonterminating 2<Pi series. With the help 
of a highly nontrivial summation formula, proved by Rahman in the appendix of jH] (see 
[7j for extensions), this eigenfunction was expressed as one of the explicit g07-solutions 
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of the Askey- Wilson second order difference operator from [3]. This eigenfunction was 
called the Askey- Wilson function in [0], since it is a meromorphic continuation of the 
Askey- Wilson polynomial in its degree. In this paper we start by reproving this result, 
now using an explicit realization of the principal series representation of U q (sl(2,C)) as 
difference operators acting on analytic functions on the complex plane. Koornwinder's 
twisted primitive element then acts as a first order difference operator, hence eigenvectors 
are easily constructed (for the positive discrete series, this was observed by Rosengren 
in [IZj)- The corresponding matrix coefficients lead to explicit integral representations 
for eigenf unctions of the Askey- Wilson second order difference operator. These matrix 
coefficients can be rewritten as the explicit g07-series representation of the Askey- Wilson 
function by a residue computation. 

We also show that for a special choice of parameter values, the Askey- Wilson function 
reduces to an elliptic analogue of the cosine kernel. This is the analogue of the classical fact 
that the Jacobi function reduces to the cosine kernel for special parameter values, see e.g. 
[TU] . We give two proofs, one proof uses an explicit expansion formula of the Askey- Wilson 
function in Askey- Wilson polynomials from [21j . the other proof uses Cherednik's Hecke 
algebra techniques from [2] and [22] ■ 

In the second part of the paper we consider the quantum group techniques for \q\ — 1. 
In this case, the approach is similar to the construction of quantum analogues of Whittaker 
vectors and Whittaker functions from ^2]- The role of g-shifted factorials, or equivalently 
g-gamma functions, is now taken over by Ruijsenaars' ^B] hyperbolic gamma function. The 
hyperbolic gamma function is directly related to Barnes' double gamma function, as well 
as to Kurokawa's double sine function, see |19j and references therein. The quantum group 
technique applied to this particular set-up leads to an eigenfunction of the Askey- Wilson 
second order difference operator for \q\ = 1, given explicitly as an Euler type integral 
involving hyperbolic gamma functions. 

The emphasis in this paper lies on exhibiting the similarities between the < q < 1 
case and the \q\ = 1 case as much as possible. Other approaches might very well lead to 
eigenfunctions for the Askey- Wilson second order difference operator for \q\ = 1 which are 
"more optimal" , in the sense that they satisfy two Askey- Wilson type difference equations 
in the geometric parameter, one with respect to base q = exp(27rir), the other with respect 
to the "modular inverted" base q — exp(27rz/r), cf. [12] for q- Whittaker functions. Such 
eigenfunctions are expected to be realized as matrix coefficients of the modular double of the 
quantum group U q (sl(2, C)) (a concept introduced by Faddeev in [3J, and are expected to 
be closely related to Ruijsenaars' ^H], 120] R-function. The R-function is an eigenfunction 
of two Askey- Wilson type second order difference operators in the geometric parameter, 
which is explicitly given as a Barnes' type integral involving hyperbolic gamma functions. 
I hope to return to these considerations in a future paper. 

Acknowledgments: It is a pleasure to dedicate this paper to Mizan Rahman. His impor- 
tant contributions to the theory of basic hypergeometric series and, more concretely, his 
kind help in the earlier stages of the research on the Askey- Wilson functions in jH], have 
played, and still play, an important role in my research on Askey- Wilson functions. 

I am supported by the Royal Netherlands Academy of Arts and Sciences (KNAW). 
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2. Generalized gamma functions 

In this section we discuss g-analogues of the gamma function for deformation parameter 
q in the regions < \q\ < 1 and \q\ = 1. 

2.1. The g-gamma function for < |g| < 1. Let r be a fixed complex number in the 
upper half plane HI. The corresponding deformation parameter q = q T = exp(27ur) has 
modulus less than one. We write q u = exp(2iriTu) for «gC. 

Let b, bj G C and n G Z> U {oo}. The g-shifted factorial is defined by 

re— 1 m 
3=0 j=l 

In g-analysis the function 

is known as the g-gamma function, see e.g. @J. For our purposes, it is more convenient to 
work with the function 



(2.i) r T (x) 



q is 



(-g^+D;g) 



Observe that T T (x) is a zero-free meromorphic function with simple poles located at — 1 + 
t -1 + 2Z< + 2r _1 Z. It furthermore satisfies the difference equation 

(2.2) T r (x + 2) = 2 cos(tt(s + l)r/2)r r (x). 



Note furthermore that for r G iM>o, i.e. < q < 1, the function r T (z) satisfies T r (x) = 
r r (^), where the bar stands for the complex conjugate. 

Observe that the above defined g-gamma type functions are not r _1 -periodic. It is 
probably for this reason that formulas in g-analysis are usually expressed in terms of g- 
shifted factorials instead of g-gamma functions. For our present purposes the expressions 
in terms of g-gamma type functions are convenient because it clarifies the similarities with 
the |g| = 1 case. 

2.2. The gamma function for |g| = 1. In this subsection we take r G IR<o, whence g = 
g T = exp(27rir) satisfies |g| = 1. As in the previous subsection, we write q u = exp(27rir-u) 
for u G C. 

It is easy to verify that the integral 

(2.3) lT {z) = i r*»(t sinh(r ^ 



2* Jo V \y sinh(?/) sinh(rj/) 
converges absolutely in the strip |Re(z)| < 1 — r _1 . For z G C in this strip we set 
(2.4) G T {z) = ex P (i lr (z)). 
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Ruijsenaars' JTHJ Sect. 3] hyperbolic gamma function G(z) = G(a + , a_; z) with a + , a_ > 
is related to G T by 

G T (z) = G(a + , a_; ia_z/2), t = —a^/a + . 

In the following proposition we recall some of Ruijsenaars' results |18[ Sect. 3] on the 
hyperbolic gamma function. 

Proposition 2.1. (i) The function G T (z) satisfies the difference equation 

G T (z + 2) = 2 cos(tt(2 + 1)t/2)G t (z). 

In particular, G T (z) admits a meromorphic continuation to the complex plane C, which we 
again denote by G T (z). 

(ii) The zeros of G T (z) are located at 1 — r _1 + 2Z>o + 2r~ 1 Z<o, and the poles of G T (z) 
are located at —1 + r _1 + 2Z< + 2r~ 1 Z> . 

(iii) G r (— z)G T (z) = 1 and G r (z) = G T -i(—rz). 

(iv) T/ie function r y T (z) has an analytic continuation to the cut plane 

C \ {(-oo, -1 + r- 1 } U [1 - r-\ oo)}, 
which we again denote by j T (z). Set r = max(l, — r _1 ) and choose e > 0, then 

2 

7TT Z 7T 

\TlAz) + — -—(T + T- 1 )\ = 0(exp((e-ir/r)\Im(z)\)), Im(z) -> ±oo 

uniformly for Re(z) in compacts o/R. 

From the explicit expression for G T (z) with |Re(z)| < 1— r _1 and the first order difference 
equation for G r , we have G T (z) = G r (z). 

As Ruijsenaars verifies in [T^ Appendix A], the hyperbolic gamma function is a quotient 
of Barnes' double gamma function, and it essentially coincides with Kurokawa's double sine 
function. 

The hyperbolic gamma function is the important building block for g-analysis with 
\q\ = 1. It was used in JB] and ^1] to construct for |g| = 1 explicit integral solutions of the 
g-Bessel difference equation and of the g-hypergeometric difference equation. Ruijsenaars' 
[TU] used hyperbolic gamma functions to construct an eigenfunction for the Askey- Wilson 
second order difference operator for |g| = 1 as an explicit Barnes' type integral. In this 
paper we construct an eigenfunction of the Askey- Wilson second order difference operator 
for |g| = 1 as an Euler type integral, using representation theory of quantum groups. 

3. The Askey- Wilson second order difference operator and quantum 

groups 

Throughout this section we require that r G C\ |Z and we write q = q T = exp(2nir) for 
the corresponding deformation parameter. The condition on r implies q ^ ±1. As usual, 
we write q u = exp(2,7[iTu) for u £ C 
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Definition 3.1. The quantum group U q is the unital associative algebra over C generated 
by K ±l , X + and X~ , subject to the relations 

KK~ l = K~ X K = 1, 

KX+ = qX + K, KX~ = q' x X~K, 

K 2 - K~ 2 



x + x~ - x~x 



q-q 1 



It is well known that U q has the structure of a Hopf-algebra, and as such it is a quantum 
deformation of the universal enveloping algebra of the simple Lie algebra sl(2, C). The 
Hopf-algebra structure does not play a significant role in the present paper, so the definition 
of the Hopf-algebra structure is omitted here. I only want to stipulate that the upcoming 
definition of Koornwinder's twisted primitive element is motivated by its transformation 
behaviour under the action of the comultiplication of U q , see ^T] for details. 

It is convenient to work with an extended version oiU q , which we define as follows. Write 
A = ©a; 6C Cx for the group algebra of the additive group (C, +). Denote End a i g (U q ) for 
the unital algebra homomorphisms : U q — > U q . There exists an algebra homomorphism 

k: A-^ End a i g (U q ), 

with k(x) = k x G End a ; 5 (W 9 ) for x G C defined by 

n x (K ±1 ) = K ±1 , K x (X ± )=q ±x X ± . 

Note that for m G Z, 

K m (X) = K m XK~ m VX G U q , 
so the automorphisms k x generalize the inner automorphisms K m ( ■ )K~ m of U q (m G Z). 
The extended algebra lA q is now defined as follows. 

Definition 3.2. The unital, associative algebra U q is the vector space A®U q with multi- 
plication defined by 

(x®X)(y®Y) = (x~T^)®K^ y (X)Y } Vi,2/eC, VX,YeU q . 

The unit element is ® 1 . 

Observe that A and U q embed as algebras in U q by the formulas 

a i-> a® 1, Inflgl 

for a G A and X G U q . We will use these canonical embeddings to identify the algebras U q 
and A with their images in U q . The commutation relations between A and U q within U q 
then become 

xX = x®X = K X (X) x, x G C, X G U q . 
The quantum Casimir element, defined by 

(q-q 1 ) i 
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is an algebraic generator of the center Z(U q ) oiU q . Note that fl is also in the center Z(U q ) 
of the extended algebra U q . 

We now consider the following explicit realization of U q . Let Ai be the space of mero- 
morphic functions on the complex plane C. For any A G C, the assignment 

±z ( q~^~ iX f(z =F 1) — q^ +iX f(z ± 1) 



{* x (X ± )f)(z) = q i 



q 1 -q 



{n x (K ±1 )f)(z) = f(z±l), 

(n x (x)f)(z) = f(z + x), xeC, 

uniquely extends to a representation of U q on A4. The quantum Casimir element fl acts 
as 

(3.2) 7r A (n) = ( ^r g C ) m - 

Observe furthermore that TT\(m) = ir\(K m ) for all m£Z. 

In the present paper the quantum group input to the theory of g-special functions is 
based on the explicit connection between the radial part of the quantum Casimir element 
fl and the second order Askey- Wilson difference operator. Here the Askey-Wilson second 
order difference operator T> = T> a,b ' c ' d , depending on four parameters (a, b, c, d) called the 
Askey-Wilson parameters, is defined by 

(3.3) (Vf)(x) = A(x)(f(x + 2) - f(x)) + A(-x)(f(x - 2) - /(*)), 

with A(x) = A(x; a, b, c, d) the explicit function 

(1 - q a+x ){\ - q b+x ){\ - q c+x ){l - q d+x ) 
y x ) ~ (1 -g 2 -)(l -g 2 + 2 -) ' 

cf. PJ. The radial part of fl is computed with respect to elements Y p — /x Q (p) 1 G U q for 
a, p G C, where Y p is Koornwinder's JT] twisted primitive element, 

(3.4) Y p = q^X + K - q^X'K + ( q P + qP \ ( K 2 - 1) 

\q 1 -q J 

and n a (p) is the constant 

'q p (l-q a )+q- p (l-q- a ) 



(3.5) fi a (p) 



q-q 1 



Consider the five dimensional space 

(3.6) U\ = span c {A+, X", K, K~\ 1} C U q . 

The radial part computation of fl leads to the following result. 
Proposition 3.3. Let p, a, a, (3 G C. For all x G C, 

xflK = xfl{x)K mod (Y p — p, a (p))xU q +xU q (Y a — p,p(a)) 
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with Q(x) = Q(x;a,p,/3,a) given explicitly by 

n(x) = ^ {B(x)K 2 + (C(x) + (1 - g 1 ^) 2 ) 1 + D(x)K- 2 } 



where 



(1 - g°+* )(l - g 2-a+x )(1 _ 6+x )(l _ 2-6+,) 

B(ar) = a, p, p,a) — q p 



(1 - q 2x )(l - q 2 + 2x ) 
C(x) = C(x; a, p, (3, a) = — A(x; a, b, c, d) — A(— x; a, b, c, d), 

D(x) = flfe a, p, /3, a) = ^ii^TX 1 " " ^H 1 " ^ 



(I - q- 2x )(l - q 2 - 2x ) 

Here the Askey- Wilson parameters (a, b, c, d) are related to the parameters a, (3, p, a by 

(3.7) (a,b,c,d) = (1 + p + cr, 1 — p + o~, 1 + a + p — (3 — a, 1 — a — p — (3 — a). 

Proof. The proof generalizes the radial part computation by Koornwinder [TT], where the 
case a = (3 = and x G Z is considered (see also Noumi & Mimachi ^3] and Koelink 
jH] for extensions to discrete values of a and (3). In the present set-up the computation is 
a bit more complex, and we gather more precise information on the remainder. For the 
convenience of the reader, I have included the main steps of the proof as appendix. □ 

Unless specified otherwise, we assume that the Askey- Wilson parameters (a, b, c, d) are 
related to the four parameters (a,p,f3,a) by ()3.7|) . 

Proposition 13.31 allows us to identify specific eigenfunctions of tt\(Q) with eigenf unctions 
of the second order difference operator 

(3.8) (Cf)(x) = (C a '^f)(x) = B(x)f(x + 2) + C(x)f(x) + D(x)f(x - 2). 

The second order difference operator £ is gauge equivalent to the Askey- Wilson second 
order difference equation D = T> a • ,c " ■ , since C = A o V o A" 1 with A(x) = A(ar; a, 6, c, d) 
any meromorphic function satisfying the difference equation 

sin(7r(2 — a + x)t) sin(7r(2 — b + x)r) 



l-g c+x )(l -g d+x ) 



(l-g2-a+x)( 1 _q r 2-6+*) 

To make use of the above radial part computation, we first need to construct explicit 
eigenfunctions of TT\(Y p ) with eigenvalue p a (p)- As we will see in the proof of the following 
proposition, the operator tt\{Y p ) is a first order difference operator, hence eigenfunctions 
of 7T\(Y p ) admit the following simple characterization. 

Proposition 3.4. Let a,p G C. A meromorphic function f G M is an eigenfunction of 
TT\{Yp) with eigenvalue p a (p) if and only if 

+ 2) — ^ — oc — p + z)t) sin(7r(-iA + a + p + z)t) 

sin(7r(iA - p + 1 + z)r) sin(7r(iA + p + 1 + z)r) 
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Proof. A direct computation shows that tt x (Y p ) G End(.M) is the explicit first order differ- 
ence operator 

(n x (Y p )f)(z) = {((f- 1 -*-* - 1)(1 - q^+ iX+z )f(z + 2) 

_ iqP+ iX-z_ 1) ( 1 _ gP -iX+z )f{z) y 

and, more generally, 

(MY p ) - = ~^~Z~q~i {W MX ~ Z ~ 1)(1 " <f +1+lA+ *)/(^ + 2) 

The eigenvalue equation 7r A (F p )/ = p a (p)f is thus equivalent to the first order difference 
equation 

^ qa+p+i\— z \fi qOL+p-i\+z^ 

f(z + 2)= qP -i-ix- z ^i _ q p+i+ix+z-j 1 a f( z )- 

Rewriting this formula yields the desired result. □ 

4. The Askey- Wilson function for < q < 1. 
In this section we take r G iM>o, so that < q = q T = exp(27rzr) < 1. The assignment 
(4.1) (K ±1 )* = K ± \ (X±)* = -X T 

uniquely extends to a unital, anti-linear, anti-algebra involution on U q . This particular 
choice of *-structure corresponds classically to choosing the real form su(l, 1) of st(2, C). 

The elements K m (m € Z), the quantum Casimir element Q (see (|3.1|) ). and the special 
family Y p (p G M.) of Koornwinder's twisted primitive elements (|3.4|) are *-selfadjoint 
elements in The eigenvalue p a (p) (see (|3.5|0 is real for a,p G R. We consider now an 
explicit *-unitary pairing for the representation n x . 

Lemma 4.1. Let AgK. Suppose that f,g £ A4 are r _1 -periodic and analytic on the strip 
{z G C| |.Re(*)| < 1}. T/ien 

(7T X (X)f,g) = (f,7T X (X*)g), WXeUl 

with the pairing (-, •) defined by 

(f,g) ■■= / f(vMg(y/T)dy. 
Jo 

Proof. This is an easy verification for the basis elements 1, K , X of U^, using Cauchy's 
Theorem to shift contours. □ 

Remark 4.2. This lemma can be applied recursively. Let /, g G M. be r _1 -periodic and 
analytic on the strip {z G C | |Re(z)| < k} with k G Z>o. For any X = X1X2 • • -X m G W g 
with m < k and Xj G ZY*, 



(n x (X)f,g) = (f,n x (X*)g). 
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In particular, the subspace T -i of entire, r _1 -periodic functions is an *-unitary n\- 
invariant subspace of M. with respect to the pairing (■,■). This subspace is the algebraic 
version of the principal series representation of the *-algebra (U q , *). 

To combine this lemma with the radial part computation of the quantum Casimir element 
(see Proposition 13.3)1 . we need to construct meromorphic r _1 -periodic eigenf unctions of 
7T\(Y p ) (p G R) which are analytic in a large enough strip around the imaginary axis. We 
claim that the meromorphic function 

(4.2) h (z) = f x (z; a, p) := M-l -£ + «» + /.- iA + z)T^(-j + p-i\-z) 

r 2T (i-i + a + p + u- z )r 2T (-i + ^+a + z) 

meets these criteria for special values of the parameters. By Proposition 13.41 and the 
difference equation for T T it follows that fx(z;a,p) is an eigenfunction of vr A (Y" p ) with 
eigenvalue p a (p) for a, p, A G C. Writing f\(z; a, p) in terms of g-shifted factorials leads to 
the expression 

( 2+a+p+iX-z 1+p+iX+z. 2\ 

/ A n\ C I \ /~i W >" ifVl _OZ 

for some nonzero constant C (independent of z), hence f\(z) is r _1 -periodic if and only if 
a G 2Z. Furthermore, observe that the poles of f\(z) are located at 

(4.4) iA-a-p + 2Z<o + Zr~ 1 , -iA + p + 1 + 2Z> + Zr" 1 , 

so f\{z) is analytic on the strip {z G C | |Re(z)| < p} when a > 0, p > and A G R. 

Definition 4.3. Let A G R, G 2Z >0 and p, a G R> 3 . For x G C wift |i?e(x)| < 2 u/e 
de/zne A (x) = A (x; a, p, /3, a) by 

4>x{x) := (TT X (xK)fx(-;/3,a)Jx(-;a,p)). 

Note that the matrix coefficient <j)\(x) is given explicitly by 

ml r*(-± + /3 + *-i\ + x + *)r* r (-l-±+<r-i\-x-iL) 



0a(z) 

(4.5) 



r 2T (-±+p + a + i\-x-z)r 2T (i-± + <j + i\ + x + z) 
r 2r (-1 + ^ + a + p + a - f )r 2r + p + zA + f ) 



r 2T (i + i + a + p - 1\ + y - )r 2r (i + p - zA - f ) 

and that 0a (#) is analytic on the strip {x G C | |Re(x)| < 2}. The quantum group inter- 
pretation of this explicit integral leads to the following result. 

Theorem 4.4. Let A G R, a, (3 G 2Z>o and p,a G R>3. T/ie matrix coefficient <px{x) = 
(f)x(x;a,p,P,a) satisfies the second order difference equation 

(Ccbx){x) = E{\)cbx{x) 

for generic x G z'R, with the eigenvalue E(X) = E(X; (3) given by 

(4.6) E{\) = -1 - q 2 ^ + q x -\q 2lX + q~ 2iX ). 
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Proof. For the duration of the proof we use the shorthand notations f(z) = f\(z; (3, a) and 
g(z) = f\(z] a, p). By the conditions on the parameters, Y p is *-selfadjoint, p a (p) is real and 
the meromorphic functions f(z) and g(z) are analytic on the strip {z G C | |Re(z)| < 3}. 
By Lemma 14.11 we thus obtain for any X G U). and x G £R, 

(n\((Y p - fi a (p))xX)f, g) = (irx(xX)f, tt x (Y p - nM)g) = 0, 

and obviously also (ir\(xX(Y a — /jLp(cr)))f, g) = 0. We conclude from Proposition 13.31 that 

(n x (xQK)f,g) = (n x (xQ(x)Kf,g) 

qP- 1 (n± ,, . q^Hl-q 1 ^) 2 



(£0 A ) (x) + — — 4>x(x). 



(9-9 1 ) (q-q x ) 2 

On the other hand, ()3.2j) implies that 

(n x (xnK)f,g)=( qlX ~ q l X ) </> A (x), 

V q-q ) 

hence (£.</>x)(x) = E{\)<f> x (x). □ 

For the comparison with the results for |g| = 1, see Section 5, it is convenient to note 
that the matrix coefficient <p x (x) can be rewritten as 

1 r 2T (± + P + a-i\ + x + I )r 2T (-l + ± + a-i\-x-^) 
r 2r (-£ + (3 + 0- + iX - x - y - )r 2r (1 - i + a + iX + x + y - ) 

x r 2r (-i-^ + ^ + p + ^-f)r 2T (-^ + p + g + f) 

r 2r (l - i + a + p - iX + y -)T 2T (-±+ p -iX-^) 

for some x-independent nonzero constant C. Formula (|4.7|) follows from ()4.5jl and the 
difference equation 

r 2r (x + — ) =exp(-^)r 2r (a;- — ). 

Corollary 4.5. Let A G R, a, /3 G 2Z>o and p, cr G M>3- Let £/ie Askey-Wilson parameters 
(a,b,c,d) be given by (13. 7|) . T/ie function F x (x) = Fx(x;a,b,c,d) defined by F x (x) = 
A(x) _1 0a(^)j wrf/i A(x) = A(x;a,b,c,d) the t~ 1 -periodic, meromorphic function 

(4.8) A W = M-l + ^ + a-x)M-l + ^ + 6-x) 

satisfies the Askey-Wilson second order difference equation 

(VF x )(x) = E(X)F x (x) 
for generic x G iM. and generic p and a. 
Proof. Note that A(x) can be rewritten as 
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for some nonzero (a;- independent) constant C . Since (3 G 2Z >0 , the gauge factor A(x) is 
r _1 -invariant. Furthermore, A(x) _1 is regular at x G ±2 + zR and x G iR under generic 
conditions on the parameters p and cr. The proof is completed by observing that A(x) 
satisfies the difference equation (J3.9jl . □ 

We end this section by extending these results to continuous parameters a and f3. Chang- 
ing integration variable and substituting the expression for FV in terms of g-shifted facto- 
rials, we can rewrite the matrix coefficient <fi\(x) (see (|4.5jl ) as 



(4.10) 6 x (x) = C- / ^ — 1 \ Joo z- 



27TZ 7 X (q^^/z, gl+/3+^-iA+x^ ) q l+P+i\ Zj q a+p+i\ / z - q 2^ ^ 

for some ^-independent nonzero constant C, where T = {z G C | \z\ = 1} is the positively 
oriented unit circle in the complex plane. To allow a,/3 to be continuous parameters, we 
need to get rid of the term z~^~ in the integrand. This can be achieved by rewriting 4>\(x) 



as 



q-¥ f (g^+'+tt-*/*, q-L+o+P-P-^/z; q 2 ) c 



d>\(x) = C 
(4.11) " { 



(ql-<x-p+P+i\ Zj q 2+a+i\+x z ^ q 2+a+p~i\ z . tfj^fa 



with C again some (different) irrelevant ^-independent nonzero constant. The integral 
formula (J4.ll)) follows from (|4.1(J|) by substitution of the identity 

Z > {q^ ^/z;q 2 ) oo =(-q^ ') > g 4 V * 

which in turn is a direct consequence of the functional equation 

6(q 2k z) = (-z)- k q- k2 6{z), keZ 

for the modified Jacobi theta function 9(z) = (qz, q/z; q 2 )^- By (14.9)1 . the eigenfunc- 
tion F\(x) = F\(x;a,b,c,d) of the Askey- Wilson second order difference equation V (see 
Corollary 14.5)1 is equal to 

f q l+p+(T—x ql-p+a-x. q2\ 

F x (x) = F x (x- a, b, c, d) := ^ ql _ a _ p+p+a J^ ql+a+p+ p +a _™ q2 ^ 

(4.12) x / ^ 7 — / - Jo ° 

V ; 2vri7 T (q^P+^z, q l +P+ iX z; q 2 ) ^ 

( q l-<x-P+P+i* Zj q 2+v+i\+x Zj q 2+a+p-i\ z . g 2^ ^ 
(q a +P+ iX /z, q a ~ iX ~ x /z, gl+/3+<T-iA+z z . q 2^ z 

up to some nonzero x-independent multiplicative constant. 
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Define for the Askey- Wilson parameters (a,b,c,d) given by (|3.7)L dual Askey- Wilson 
parameters (a, b, c, d) by 

(a,b,c, d) = (1 - P, 1 + /3 + 2(7, 1 + a + 2p, 1 - a). 

This notion of dual Askey- Wilson parameters coincides with the notion of dual parameters 
as used in e.g. 0. 

In the following theorem we express T\{x) in terms of basic hypergeometric series by 
shrinking the radius of the integration circle T to zero while picking up residues. Recall 
that the very-well-poised 807 series is defined by 

(I - uq 2k )(u,b 1 ,b 2 ,b 3 ,b4,b 5 ;q) z k 

8 W 7 (u; h, b 2 , b 3 , 6 4 , h- q,z)=}_^ - r 

^ (1 -u){q,qu/b 1 ,qu/b 2 ,qu/b 3 ,qu/h,qu/b 5 iq) k 

for \z\ < 1, see jl]. 

Theorem 4.6. Let A 6 1, «J G R>o and p,o~ G M>3. Under the parameter correspon- 
dence (|3.7j) . t/ie function F\{x) = J"a(x; a, 6, c, d) given 6?/ (j4.12j) can fee expressed in terms 
of basic hypergeometric series as 

I 2+a-d+2i\+x 2+a-d+2iX-x. 2\ 



(4.13) (g 2 -^,?^--? 2 )^ 

X s W 7 (q- 2+ ~ a+ ~ b+ ~ c+2lX - q a+x , q a ~ x , g S+2i \ q l+2l \ q Z+2l \ q 2 , q 2 -^) 
with the [irrelevant] generically nonzero, x-independent constant C given by 

( n 2+2a 2+a+2p-/3 2+a+2p+/3+2cr l+/3+2iA l+a-2iA. 2\ 
^1 1 H 1 H 1 H ) H '"/oo 



C 



( n 2 n l+a+2i\ n l+a+2p+2i\ n l+(3+2a-2iX n 3+a+2p+2a+2i\. n 2\ 



Furthermore, J-\{x) is t 1 -periodic and satisfies the Askey-Wilson difference equation 
(T>J-'xj(x) = E(X)J-'\(x) for generic x G iM. 

Proof. The expression for T\ follows by shrinking the radius of the integration contour T 
to zero while picking up residues. It is actually a special case of fU Exerc. 4.4, p. 122], in 
which one should replace the base q by q 2 and the parameters (a, 6, c, d, f, g, h, k) by 

( l-0-a-i\+x 1-p+iX a+p+iX a-iX-x 2+a+iX+x 1+p+iX l+/3+a-iX+x l-/3-2«A\ 

yq 1 q 1 q 1 q 1 q 1 q 1 q 1 q ) ■ 

We have seen that the difference equation (DjF A )(x) = £'(A)jF A (x) for generic x G iM. is 
valid under the extra assumption a,j3 G 2Z >0 . For a,/3 G !R>o this difference equation has 
been proved by Ismail and Rahman [5 using the explicit expression of T\{x) as very- well- 
poised 807 series. □ 

Using the explicit expressions of J-"\(x), the conditions on x, A and the four parame- 
ters a,p,/3,o~ can be relaxed by meromorphic continuation. The resulting function T\ is 
a meromorphic, r~ Aperiodic eigenfunction of the Askey-Wilson second order difference 
operator V = V a ^ c4 with eigenvalue E{\) = E(X; (3). 
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Remark 4.7. Some special cases of the matrix coefficients <p\ were explicitly expressed 
in terms of very-well-poised §07 series in jS] using the realization of the representation 
tt\ on the representation space / 2 (Z). In this approach the basic hypergeometric series 
manipulations are much harder, since one needs a highly nontrivial evaluation of a non- 
symmetric Poisson type kernel involving nonterminating 20i-series which is due to Rahman, 
see the appendix of [S] and jjj. 

Remark 4.8. The explicit §07 expression (j4.13|) of T\(x) and its analytic continuation 
was named the Askey- Wilson function in 0. Suslov [21] established Fourier-Bessel 
type orthogonality relations for the Askey- Wilson function. Koelink and the author jS] , [5| 
defined a generalized Fourier transform involving the Askey- Wilson function as the integral 
kernel, and established its Plancherel and inversion formula. This transform, called the 
Askey- Wilson function transform, arises as Fourier transform on the noncompact quantum 
group SU ? (1, 1) (see [8J, and may thus be seen as a natural analogue of the Jacobi function 
transform. 



5. The expansion formula and the elliptic cosine kernel 

In this section we still assume that r G iM>o, so < q = q T = exp(27iir) < 1. To keep 
contact with the conventions of the previous section, we keep working in base q 2 . 

First we recall the (normalized) Askey- Wilson polynomials [I]. The Askey- Wilson poly- 
nomials are defined by 

/„-2m Jim— 2+a+b+c+d „a+x „a—x \ 

E m (x) = E m (x] a, b, c, d) := 4 3 ' +b +c ' +d ' ; q 2 , q 2 , me Z> , 

\ H 1 H 1 H J 



with 



ai,a 2 ,a3,a4 \ _ (°i> a 2; a 3, 0,4, q) k k 



m h , , ;q,z 



E^u,i, u,^, u, 3 , u, 4 , H j k k 111 
'/ \ & 1 z J- • 



' k 

The Askey- Wilson polynomial E m (x) is a polynomial in q x + g _I of degree m, normalized 
by E m (a) = 1. They satisfy the orthogonality relations 

1 (q 2x , q~ 2x : q 2 ) 
E m (x/T)E n (x/T) -, ^—^ Jo ° - dx = 0, m^n, 

v ' ' \ 1 / I „a+x n a—x n b+x n b—x n c+x n c—x n d+x n d—x- n l\ ' 
yd -i h 1 H 1 h 1 H 1 H 1 h 1 H ' " / 00 

provided that Re (a), Re(6), Re(c), Re(d) > 0. Observe that the Askey- Wilson polynomial 
E m {x) is regular at the special choice 

(5.1) (a,b,c,d) = (0,-L 1,1-^-) 

of Askey- Wilson parameters. The above orthogonality relations also extend (by continuity) 
to the Askey- Wilson parameters (j5.1j) . leading to 

f 1 11 11 

/ E m (x/T;0,—,l,l-—)E n (x/T;0,-,l,l-—)dx = 0, m ^ n, 

J 
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hence we conclude that 

(5.2) E m (x;0, — , 1, 1 - -!-) = cos(27rmrx), m G Z> , 

zr zr 

which is the usual cosine kernel from the Fourier theory on the unit circle. In this section 
we derive an analogous result for the Askey- Wilson function. 

In the previous section we have introduced the dual Askey- Wilson parameters associated 
to a, b, c and d. They can be alternatively expressed as 

1 1 1 1 

(a, b, c, d) = (—l + -(a+b+c+d),l + ~(a+b—c—d),l + -(a—b+c—d),l + -(a—b—c+d) > ). 

Note that the special choice (|5.1|) of Askey- Wilson parameters is self dual, (a, b, c, d) = 
(a,b,c,d). For our present purposes it is convenient to use yet another normalization of 
the Askey- Wilson function, namely 

x) = x; a, b, c, d) 

( a 2+a— d+iM+x n 2+a—d+fi—x J2—a—d „2+o— d. „2\ 
\H >_y ij_ >_* ' " I oo 

^gd+b+c+fj, g2—d+fi g2—d+x g2—d—x- q2^ 

x % W 7 {q- 2+ ~ a+l+ ~ c+ ^- q a+x , q a ~ x , q^, q b+ ^ q^- q\ g 2 "^), |g 2 ~^| < 1. 

For fixed A, the eigenfunction of the Askey-Wilson second order difference operator 

D is a constant multiple of <£ + (2zA, •). The present normalization of the Askey-Wilson 
function is convenient due to the properties 

(5.3) ^ + (/ i > x \ a i b, c, d) = £ + (x, \i\ a, b,c,d) 

and (E + (— a, — a) = 1. The property ()5.3j) is called duality and can be proved using a 
transformation formula for very-well-poised §07 series, see [5] for details. Furthermore, 

(5.4) (£ + (a + 2m, x) = E m (x), m e Z> , 

see e.g. (3.5)], thus the Askey-Wilson function (£ + (/i, x) provides a natural meromorphic 
continuation of the Askey-Wilson polynomial in its degree. 

The meromorphic continuation of £ + (/i, x) in /i and x can be established by the integral 
representation of the Askey-Wilson function (see the previous section), or by the expression 
of the Askey-Wilson function as a sum of two balanced 403's (see e.g. (3.3)]). For our 
present purposes, it is most convenient to consider the meromorphic continuation via the 
expansion formula of the Askey-Wilson function in Askey-Wilson polynomials, given by 

( n 2— a— d Jl+a—d „b+c „2+b—d J2+c—d. „2\ 

\q i q i q i q i q '"/oo 



£+(/!, X) 



( n 2+a+b+c— d n 2—d+x n 2—d—x n 2—d+fi n 2—d—fi. n 2\ 
\H i H i H i H i H i y J j 



(5.5) x E m (x; a, b, c, 2 — d)E m ([i] a, b, c, 2 — d ) 

m=0 

1 1 n 2m+a+b+c— d\ I n a+b+c—d „a+b n a+c. n 2\ 

x y 1 1 >J± )m f (l-q-ri)ro m? 

(l_ g a+6+ C -d)( ? 2 )g 2+6-d ;? 2 + c-d ;g 2) * ' 



ASKEY-WILSON FUNCTIONS AND QUANTUM GROUPS. 15 

see |2*T1 Thm. 4.2]. The sum converges absolutely and uniformly on compacta of (/J,,x) G 
C x C due to the Gaussian q m . The expansion formula ()5.5|) shows that the Askey- Wilson 
function <£ + (/i, x) is well defined and regular at the special choice (|5.1|) of Askey- Wilson 
parameters. In fact, for this special choice of parameters, the Askey- Wilson function can 
be expressed in terms of the (renormalized) Jacobi theta function 

#(x) = {-q 1+ *,-q 1 -*rf) oo 

as follows. 

Proposition 5.1. We have the identity 

Proof. To simplify notations, we write 

(5.7) (£+(/i, x) = (E + (n, x- 0, i-, 1, 1 - 

for the duration of the proof. 

We substitute the special choice (|5.1|) of Askey- Wilson parameters in the expansion 
formula for (£ + (/i, x). By simple (/-series manipulations and by (|5.2|) . we obtain the explicit 
formula 

x) = — 1 + 2 > cos(27rmr/i) cos(27rmrx) g 

UWj^WW V m=l 

Using the well known Jacobi triple product identity 

oo oo 

J2 q m> + m X = i + 2 J2 cos(2nmrx) q m * = (g 2 ; q 2 ) ^ 0(x), 

m=— oo m=l 

and the elementary identity 

cos(27rmr/i) cos(27rmrx) = - (cos(27rmr( / u + x)) + cos(27rmr(/i — x))) , 
we deduce that 

+ , , _ (-g.-g.-g^-gWL f ti{fj, + x) + 0{(j.-x) 

olM ' j " 2 V ^)^(x) 

Simplifying the multiplicative constant yields the desired result. □ 

Remark 5.2. The Jacobi theta function $(x) is the natural r _1 -periodic analogue of the 
Gaussian q~ x , since 

oo 

£ g-^^v^^g 2 )^*) 

fc=— oo 

by the Jacobi triple product identity and the Jacobi inversion formula. In fact, in j2Ij 
and [22] it is shown that the function (q 2 ~ d+x , q 2 - d ~ x ^ g 2 )^, which reduces to $(x) for the 
Askey- Wilson parameters (|5.1|) . plays the role of the Gaussian in the Askey- Wilson theory. 
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If one replaces the theta functions by Gaussians in the right hand side of (|5.6|h then we 
obtain up to a multiplicative constant 

' -(/i+a:) 2 _|_ n -(n-x) 2 \ 1 

j = i(g"^ + 9*") = cos(47rr^), 

which is essentially the classical cosine kernel. Thus the right hand side of ()5.6|) is an 
elliptic analogue of the cosine kernel. 

Remark 5.3. Using the quasi-periodicity 

(5.8) •&{% + 2) = ^-^(a;) 

of the Jacobi theta function, we obtain as a consequence of 1)5. 6 J) . 

(E + (2m, x; 0, — ,1,1 ) = cos(27rmrx), m G Z> , 

It It 

which is in accordance with ()5.2|) and (|5.4|) . 

Remark 5.4. The orthogonality relations for the Askey- Wilson polynomials with Askey- 
Wilson parameters (|5.1j) are equivalent to the L 2 -theory of the classical Fourier tansform on 
the unit circle. On the other hand, the L 2 -theory of the Askey- Wilson function transform, 
see jH], does not reduce to the L 2 -theory of the classical Fourier theory on the real line 
for the Askey- Wilson parameters (|5.1j) . Instead one obtains a Fourier type transform with 
integral kernel given by the elliptic cosine function ([5.6)1 . For the corresponding L 2 theory, 
the transform is defined on a weighted L 2 -space consisting of functions that are supported 
on a finite closed interval and an infinite, unbounded sequence of discrete mass points. 
This transform, as well as the general Askey- Wilson function transform, still has many 
properties in common with the classical Fourier transform on the real line, see e.g. [Oj, [21] 
and [22] 

Cherednik's j2j Hecke algebra approach to g-special functions leads to a direct proof 
that the right hand side of the expansion formula ()5.5|) is an eigenfunction of the Askey- 
Wilson second order difference operator T>, see [22J. The expansion formula ()5.5|) may thus 
be seen as the explicit link between Cherednik's approach and Ismail's and Rahman's [S] 
construction of eigenfunctions of T> in terms of very-well-poised g</>7 series. We end this 
section by sketching a proof of Proposition 15. II using Cherednik's Hecke algebra approach. 

The affine Hecke algebra techniques for Askey- Wilson polynomials are developed in full 
detail in ^H] , and for Askey- Wilson functions in [22] • We first recall one of the main results 
from [22], specialized to the present rank one situation. 

Define two difference-reflection operators by 



(5.9) 



(T c /f){x) = -q-^ +d f{x) + (1 t X J { l-S X) (/( 2 - *) - /(*))> 



(T?> b f)(x) = - q ^ b f(x) + (1 q ^2 X) qb+X) (/("*) - /(*))■ 
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The connection with affine Hecke algebras follows from the fact that T = Tq ' d and T\ = T® ,b 
satisfy Hecke type quadratic relations. These relations imply that the operators T and Ti 
are invertible. Consider the (invertible) operator 

(5.10) Y = Y a ' b ' c ' d := T? h o T c ' d . 

Remark 5.5. The operator Y + Y^ 1 , acting on even functions, is essentially the Askey- 
Wilson second order difference operator V, see e.g. [TBI Prop. 5.8]. 

Theorem 5.17 in (221 states that for generic Askey- Wilson parameters (a,b,c,d), there 
exists a unique meromorphic function <£(•,•) = £(•, •; a, b, c, d) on C x C satisfying the 
following six conditions: 

1. (£(//, x) is r _1 -periodic in /x and x, 

2. (fijX) i— >• (q 2 ~ d+ L\q 2 ~ d ~^,q 2 ~ d+x ,q 2 ~ d ~ x ]q 2 ) ^^[i^x) is analytic, 

3. For fixed generic /i e C, (£(//, •) is an eigenfunction of Y a,b,c ' d with eigenvalue g a_/i , 

4. For fixed generic x G C, <£(-,zc) is an eigenfunction of y a ' b > c ' d with eigenvalue q a ~ x , 

5. (T^(/i,-))(x) = —q a+b (Tf'*€(-, 

6. (£(— a, —a) = 1. 

The existence of a kernel (£ satisfying the above six conditions is proved by explicitly 
constructing (E as series expansion in nonsymmetric analogues of the Askey- Wilson poly- 
nomials, see |2~2"l (6.6)]. This expansion formula for <£ is very similar to the expansion 
formula ()5.5|) of the Askey- Wilson function (£ + in Askey- Wilson polynomials. In fact, a 
comparison of the formulas leads to the explicit link 

(5-11) <£ + (/i,o;) = (C+<£(/v))(s), := (l +T?> b ), 

see |22l Thm. 6.20]. These results allow us to study the Askey- Wilson function (E + using 
the characterizing conditions 1—6 for the underlying kernel (£, instead of focussing on the 
explicit expression for 

The kernel (£(//, x) is regular at the Askey- Wilson parameters (|5.1j) . The resulting kernel 

Co(//,a;) = <£(//, z; 0,-^,1,1--!-), 
zr zr 

is the unique meromorphic kernel satisfying the six conditions 1—6 for the special Askey- 
Wilson parameters (J5.1)) . Observe that the operators T ,Ti and Y for the special Askey- 
Wilson parameters (|5.1j) reduce to 

(T f)(x) = f(2 - x), (Tx/)(x) = /(-*), (Yf)(x) = f(2 + x), 
hence (Eq is the unique meromorphic kernel satisfying the six conditions 
V. <£ (fi,x) is r _1 -periodic in \i and x, 
2'. (/i,x) i — > i?(/i)i?(x)(£o(A*, is analytic, 
3'. e (^,x + 2) = g-"C (//,a;), 
4'. £0(^ + 2,^) = 9" x ^o(/i,a;), 
5'. (£ (/i, -x) = (£ (-/i,a;), 
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6'. <* (0,0) = 1. 

We conclude that 

£o(/,,o:) = H,-g;g)^^ y ^, 

since the right hand side satisfies 1'— 6' due to the quasi-periodicity ()5.8|) of -3(x). Thus 
(£q(jU, x) is an elliptic analogue of the exponential kernel exp(— Arrirpx), cf. Remark 15.21 
Using the notation (|5.7jl . we conclude that 

£o(v,x) = (C+^e (/i, •))(») = gC^oO^aO + g o(/i, -ar)), 
which is the desired formula (|5.6jl . 

6. The Askey- Wilson function for |g| = 1. 

In this section we take — | < r < 0, so that q = q T = exp(2irir) has modulus one and 
q ^ ±1. The assignment 

(6.1) (if* 1 )* = IT 41 , (X ± )' t = -X ± 

uniquely extends to a unital, anti-linear, anti-algebra involution on U q . This particu- 
lar choice of ^-structure corresponds to the real form st(2, R) of st(2, C). The ^-unitary 
sesquilinear form for the representations 7r A (A G R) of W g (see Section 3) is 

/ioo 
f{z)g(z) dz, 
■ioo 

where / and g are meromorphic functions which are regular on a large enough strip around 
the imaginary axes and decay sufficiently fast at ±zoo. Koornwinder's twisted primitive 
element iY p 6 U q is *-selfadjoint for p e R. Thus in principle we are all set to extend the 
construction of eigenfunctions of the Askey- Wilson second order difference operator T> to 
the \q\ — 1 case by simply replacing the role of the g-gamma function IV by G2 T - We need 
to be careful though due to the following differences with the < q < 1 case: 

a. The analogue of the explicit eigenfunction of iY p (cf. fl4.2j) ). given now as quotient 
of hyperbolic gamma functions G^t, has more singularities. 

b. No r^ 1 -periodicity conditions have to be imposed. Consequently, the parameters 
a and j3 do not need to be discretized for the \q\ = 1 case. 

c. We have to take the decay rates at ±ioo of integrands into account. 

One needs to be careful with the decay rate (see c) in reproving the crucial Theorem 14.41 
because acting by 7Ta(A ± ) worsens the asymptotics at ±ioo (the factor q z = exp(2niTz) 
is 0(exp(—2iTTlm(z)) as Im(^) — > oo). The decay rate can be improved by considering 
different eigenfunctions of iY p , but then the location of the singularities turns out to cause 
problems. 

To get around these problems, we generalize the techniques of Section 4 to a nonunitary 
set-up. More concretely, we replace (•,•)' by a bilinear form, given as a contour integral 
over a certain deformation of zR. With such a bilinear form, the singularity problems and 
the asymptotic problems can be resolved simultaneously. 
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The proper replacement of the involution * is the unique unital, linear, anti-algebra 
involution o on U q satisfying 

{K ±l )° = K^\ {X^ = ~X ± . 

For / G Ai we write S(f) C C for the singular set of /. 

Definition 6.1. / G M. is said to have (exponential) growth rate e G R at ±ioo when the 
following two conditions are satisfied: 

1. For some compact subset Kf C R, 

S(f) C {z G C | /m(z) G K/}. 

2. T/ie function f satisfies 

\f(x + iy)\ = 0(exp(e\y\)), y -> ±oo, 

uniformly for x in compacts o/R. 

We call a contour C a deformation o/iR when C intersects the line 

/ c = {2; G C | Im(z) = c} 

in exactly one point z c (C) for all c G R, and z c (C) = zc for |c| >> 0. For G Z> we define 
the strip of radius k around C by 

|J{^ G / c I \z-z c (C)\ < k}. 



For k = 0, the strip of radius zero around C is the contour C itself. 

Lemma 6.2. Suppose that f,g E Ai have growth rates ej and e g at ±ioo, respectively. 
Suppose furthermore that e/ + e 9 < 2ttt and that f and g are analytic on the strip of radius 
one around a given (oriented) deformation C o/iR. Then 

(n x (X)f,g) c = (f,^ x (X°)g) c , VIgWJ, 

where 

(f>9) c = / f( z )d(z)dz. 



Proof. The proof follows by an elementary application of Cauchy's Theorem. □ 

Remark 6.3. The condition on the growth rates in Lemma lfi~21 mav be weakened to tf + t g < 
when X is taken from the subspace span c {l, K^ 1 , K} of UL 

Koornwinder's twisted primitive element Y p (see (13.4(0 is not o-invariant, 
Y ; = - q kir x X+ + q-^K-'X- + (^T3y) (K- 2 - 1). 

On the other hand, a direct computation shows that n_\{Y°) is still a first order difference 
operator. This leads to the following analogue of Proposition 13.41 
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Proposition 6.4. Let a,p G C. A meromorphic function f G A4 is an eigenfunction of 
7v_\(Y°) with eigenvalue p a (p) if and only if 

j., ^ _ sm(n(iX + a + p - z)r) sin(7r(-zA + a + p + z)r) 
sin(7r(zA + p — 1 + z)t) sin(7r(— iX + p + 1 — z)t) 

We define two meromorphic functions by 

G 2 r{jz- I + P + a - iX + z)G 2r (± + a - iX - z) 



(6.2) 



g x (z; (3, a) 
h x (z;a,p) 



G 2r (-^ + l + (3 + a + i\- z)G 2r (-± + o + iX + z) 
G 2r (-^ - l + a + p + i\- z)G 2t (-± + p + i\ + z) 



G 2t (-^: + l + a + p-i\ + z)G 2t (-^ + p - iX - z) 
The difference equation for G 2r , Proposition 13.41 and Proposition 16.41 imply 

(6.3) 7r x (Y a )gx(-; f3, a) = p,p(a)gx(-; (3, a), 7r_ A (F p °)/i A (-; a, p) = p a (p)h x (-; a, p). 

We want to construct now eigenfunctions of the gauged Askey- Wilson second order differ- 
ence operator L which are of the form 

ipx(x; a, p, (3, a) = (nx(xK)g x (-; (3, a), h x (-; a, p)) c 

(6.4) f 

= J g x (l + x + z; (3, a)hx(z; a, p)dz 

for a suitable deformation C of zR. To make sense of this integral, we need to take the 
singularities and the asymptotic behaviour at ±ioo of the integrand into account. The 
singularities can be located using the precise information on the zeros and poles of the 
hyperbolic gamma function G T , see Proposition 12. llf ii) . It follows that the singularities of 
z i — > gx{l + x + z; (3, a) are contained in the union of the four half lines 

-l-P-<r + i\-x + R< , -l+P + a + i\-x + R< , 
— a — iX — x + R>o, a — iX — x + R>o, 

and the singularities of z i— > hx(z; a, p) are contained in the union of the four half lines 

- 1 + - -p -iX + R< , -l + p-zA + R< , 

r 

-a-p + iA + R> , -- + a + p + iX + R >0 . 

T 

We call the above eight half lines the singular half lines with respect to the given, fixed 
parameters r,x, A, a, p, P,a. Each singular half line is contained in some horizontal line 
l c for some c G R. For generic parameters a,p,P,o~, the eight singular half lines lie on 
different horizontal lines. Under these generic assumptions, there exists a deformation C x 
of zR which seperates the four singular half lines with real part tending to — oo from the 
four singular half lines with real part tending to oo. We take such contour C = C x in the 
definition ()6.4|) of ij) x - The resulting function ipx(x) is well defined and independent of the 
particular choice of the deformation C x of iR, since g x (-', P, cr) (respectively hx(-] a, p)) has 
growth rate 7r((l — 2Im(A))r — l) (respectively tt(1 + 2Im(A))r) at ±ioo. This follows from 
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the asymptotic behaviour of the hyperbolic gamma function G T , see Proposition 12. ll fiv) . 
The resulting function ip\{x) is analytic in x. 

Theorem 6.5. For generic parameters a,p,j3,a, the function = a, p, (3, a) de- 
fined by 



(6.5) 



ipx(x) =/ gx(l + x + z;P,a)hx(z;a,p)dz 
>c x 

G 2T (± + f3 + a-iX + x + z)G 2t (± -1 + a-iX-x-z) 
"c G 2t + P + o- + iX - x - z) G 2r (-± + l + a + iX + x + z) 
< G 2t {-± -1 + a + p + tX- z)G 2t {-± + P + iX + z ) d 
G 2t + 1 + a + p - iX + z)G 2t (-i + p - iX - z) 

is an eigenfunction of the gauged Askey- Wilson difference operator C a,p, ^ ,cr with eigenvalue 
E(X;(3). 

Proof. We adjust the proof of Theorem 14.41 to the present set-up. We simplify notations 
by writing g(z) = g\(z; (3, a) and h(z) = h x (z; a, p). Choose the deformation C x of iR such 
that g and h are analytic on the strip of radius > 4 around C x . 

Observe that n x (xQK)g has the same growth rate at ±ioo as g. By (|3.2|) and the 
definition of tp\(x), we conclude that the integral (ir\(xT!,K)g, h) c converges absolutely 
and equals 



q-q- 1 



tpx(x). 



On the other hand, the radial part computation of f2 with respect to Koornwinder's twisted 
primitive element yields 

(6.6) xVLK = m(x)K + ((Y p - p a (p))K^)XxK + xZ(Y a - p, 3 {a)) 

for certain elements X, Z gW, 1 , cf. Proposition 13.31 Substituting this algebraic identity in 
(n x (xQK)g, h) c and using 7r x (Y a )g = p^{a)g, we have 

(6.7) (7i x (xQK)g,h) Cx = (7i x (xQ(x)K)g,h) Cx + (tt x ((Y p - p^K^Xx^g, h)^, 

provided that both integrals on the right hand side of (|6.7j) converge absolutely. 

For the second term on the right hand side of (j6.7jl . observe that the sum of the growth 
rates of g and h at ±zoo equals (2r — 1)tt. Furthermore, 

(2r - l)n < Aitr 



since — \ < r < 0, and 



(Y p -p a (p))K- 1 eu l q . 
Hence the integral 

(tt a ((F p - p a {p))K- l XxK)g, h) Cx 
converges absolutely, and Lemma [6.21 shows that this integral equals 

(n x (XxK)g,TV^(K)n^(Y p - p a (p))h) c = 0, 
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where the last equality follows from the eigenvalue equation ir_x{Y°)h = p a (p)h. 

For the first term on the right hand side of (|6.7j) . observe that n\(xQ(x)K)g has the same 
growth rate at ±ioo as g. Hence the integral (ir\(x£l(x)K)g, h) c converges absolutely. By 
the definition of the gauged Askey- Wilson second order difference operator C, this integral 
equals 



(?-? ') (q-q i y 

Combining the results, we conclude that 

(■K X (xttK)g,h) Cx = (ir x {xYL(x)K)g,h) Cx , 

and for both sides we have obtained an explicit expression in terms of ip\- The resulting 
identity for ip\ yields the desired difference equation £ip\ = E(\)ip\. □ 

Remark 6.6. The eigenfunction ipx of C as given by the integral (J6.5|) looks very similar to 
the eigenfunction <p\ of £ (for the case r e iR>o) as given by the integral ()4.7|) . since the 
integrand of <fi\ is essentially the integrand of ip\ with the hyperbolic gamma functions £t2t 
replaced by g-gamma functions IV- Note though that the integration cycles are different. 

We can reformulate the difference equation for ip\ in terms of the Askey- Wilson second 
order difference operator T> = X> a,b,c,d (see ()3.3j0 using an appropriate gauge factor, cf. 
Corollary 14.51 Using the parameter correspondence (|3.7|) . we can for instance choose the 
gauge factor 5(x) = 5(x; a, p, /3, a) by 

( 6 .8) *(,) = + £ + \±k±^zA . 

G 2T (l + ±-c-x)G 2T (l + ±-d-x) 

Corollary 6.7. For generic parameters a, p, f3, a, the function H\ = H\(-; a, b, c,d) G A4 
defined by H\(x) = 5(x)~ 1 tp\(x) satisfies the Askey-Wilson second order difference equation 

(VH x )(x) = E(X)H x (x), 

where the Askey-Wilson parameters (a,b,c,d) are given by (|3.7|) . 

7. Appendix 

In this appendix we sketch a proof of Proposition 13.31 following closely the arguments 
of Koornwinder ^T] for the special case a = (3 = 0. 
Fix x e C and write X = X' for X, X' e U q if 

x-x'e(Y p - fi a {p)) xu\ + x wj(n - pM). 

In order to simplify notations, I use the notations p = p a {p) and v = pp(a) for the duration 
of the proof. To reduce xQK = x KQ in the desired form, we need to concentrate on the 
part x KX + X~ . Using 



x 



KX + X~ =q^ +x (q^X + K - q-^X~K)xX~ 

+ q-^ +x X~x{q^X~K - q^X + K) + q 2x xKX~X' 
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and the commutation relation between X~ and X + in U q , we obtain 



K- 1 - K 3 
q-q- 1 



(1 - q 2x ) xKX + X- =q^ x ( 9 " ' + Q " (K 2 - 1) + fi) x X~ 

\q-q J 

- q^ +x X~x (^-±ll(K 2 - 1) + i^j + q 2x ( 

hence we need to focus now only on the reduction of xX~ and xK 2 X~ . Using 

xX~ =q^~ x {q-^X'K - q^X+KjK^x 

+ q^- 2x K~ x x (q^X + K - q~^X~K) + q x ~ 2x K~ x x X~ K , 



we obtain 



(1 - q 2 ~ 2x ) xX~ = - q^- x ( ^±^ (K 2 - 1) + ^ K~ x x 

+ q^ 2x K~ x x (^-±£^(K 2 - 1) + ^ , 

and a similar computation yields 

(1 _ g-2-2x) xK 2 X~ EE - g-i"* ^ ^±i^ (X 2 - 1) + ^ i^X 



+ q-l- 2x Kx - (if - 1) + i/ . 

V ?-? / 

Collecting all these results we arrive for generic i G C at a formula of the form 

£Qif ee x (B(x)K 3 + C(x)if + Z^if- 1 ) 

for explicit rational expressions -B(rc), C(x) and D(x). Using 

g p + g~ p - (g - g" 1 )^ = g p+a + g- p - a , 
<f + q~ a - (q - <T V = <f +/3 + q~ a ~ p , 

the rational functions -B(x), C(x) and are explicitly given by 

(q 1 2 +x (qf + g"^) - gl+ 2a; (g CT + g~ CT )) (g~f- 22: (g CT + g" CT ) - q~l" x (qP + g"")) 



(g-g- 1 ) 2 (l-g 2 -)(l-g- 2 - 2 -) 



2r — 1 

g g 



(!-r 1 )(i-! h ) (g-?- 1 ) 2 
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C(X) = 



+ 
+ 

and 



(q-q- 1 ) 2 

(q^ +x (qP + q-P) - q^+ 2x (q cr + q-°)) (q-l- x (qP +a + q~P~ a ) - q-\~ 2x ((f+P + q- ff ~P)) 

{q-q- 1 ) 2 ^ -q 2x ){l - q~ 2 ~ 2x ) 
(q-\+ 2x {q lT ^ + q-^-l 3 ) - q\+ x (qP+ a + <T^ a )) (g§- 2:c (g a + q- a ) - qh~ x {qP + q~P)) 

(g-g- 1 ) 2 (l-g 2 -)(l-g 2 - 2 -) ' 



(g-g^Xi-g 2 *) (g-g- 1 ) 2 

{ql- x (qp +a + g-^ Q ) - q%~ 2 *{q°+P + 



(g-g-^^l-g^Xl-g 2 - 2 *) 
By a tedious computation, it can now be proven that 

(g-g x ) 2 

= f g ' (-A(s) - A^" 1 ) + (1 - g 1 ^) 2 ), 

with the parameters a, 6, c, cf as in ()3.7jl . 
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